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Abstract 

The evolution of the gauge coupling is modelled using a scalar degree of free- 
dom whose homogeneous and inhomogeneous modes evolve from an ordinary 
inflationary phase to the subsequent epochs. Depending upon various param- 
eters (scalar mass, curvature scale at the end of inflation and at the onset of 
the radiation epoch), the two-point function of the magnetic inhomogeneities 
grows during the de Sitter stage of expansion and, consequently, large scale 
magnetic fields are generated. The requirements coming from inflationary 
magnetogenesis are examined together with the theoretical constraints stem- 
ming from the explicit model of the evolution of the gauge coupling. Galactic 
magnetogenesis is possible for a wide range of parameters. Inter-galactic and 
inter-cluster magnetogenesis are also discussed. 
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I. INTRODUCTION 



The gauge coupling has to be (almost) constant by the time when the Universe is ap- 
proximately old of one second, namely by the time of big-bang nucleosynthesis (BBN). The 
abundances of light elements are very sensitive to any departure from the standard cos- 
mological model. Hence, fluctuations in the baryon to photon ratio ffl, matter-antimatter 
domains M, anisotropics in the expansion of the four space-time dimensions [H], can all be 
successfully constrained by demanding that the abundances of the light elements are cor- 
rectly reproduced. Following the same logic the variation in the gauge couplings can also be 
constrained from BBN [|J. 

The variation of the gauge coupling may also produce effects upon the epoch at which the 
baryon asymmetry was formed |5]|| namely the epoch of the electroweak phase transition 
(EWPT) occurring when the Universe was approximately 10~ n sec old. The influence 
of the variation of the gauge couplings on the generation of the baryon asymmetry of the 
Universe (BAU) is not direct but it can occur in a specific model. Suppose, for instance, 
that the gauge coupling evolution is modelled using a scalar field. If the field decays prior 
to the EWPT the BAU generated after the EWPT will be preserved. However, if the field 
decays between the EWPT and BBN the generated entropy associated with the decay of the 
field will have to be constrained. Too much entropy could dilute the BAU and, therefore, a 
constraint on the produced entropy will translate into a constraint on the scalar mass. 

In || it was argued that the evolution of the gauge couplings may have interesting 
effects associated with the generation of magnetic fields. In the present investigation the 
compatibility of the evolution of the gauge coupling with an ordinary (i.e de Sitter or quasi-de 
Sitter) stage of inflation will be analyzed. In this sense the present discussion is an extension 
and a completion of the results reported in || where the compatibility of the evolution of 
the gauge coupling with a de Sitter stage of inflation has not been addressed. 

The evolution of the Abelian coupling will be mainly investigated. The non-Abelian 
couplings will be fixed. This choice is motivated by the remark that the effects associated 
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with large scale gauge fluctuations are related more with the magnetic component of Abelian 
fields [[J. Furthermore, the necessity of not erasing the produced BAU (by excessive entropy 
production) leads to exclude a variability of the gauge coupling during the EWPT taking 
place (roughly) at 100 GeV [[/J. To relax these to assumptions is possible but it is beyond 
the aim of the present discussion. 

Suppose that during a de Sitter stage of expansion the coupling constant of an Abelian 
gauge field evolves in time. Thus the kinetic term of the gauge field can be written (in four 
space-time dimensions) as 

S em = ~\J d 4 xV^Gf^)F a pF a P, (1.1) 

where G is the determinant of the space-time metric, is a scalar field (which can depend 
upon space and time) and g(<fi) = /(0) _1//2 is the coupling 0. This type of vertex is typical 
of scalar-tensor theories of gravity [K| and of the low energy string effective action ||TT] . 
Early suggestions that the Abelian gauge coupling may change over cosmological times were 
originally made by Dirac [|I^] (and subsequently discussed in |T3] and in [|rj]]) mainly in the 
framework of ordinary electromagnetism. 

The dynamics of the field <fi will be described by the action of a minimally coupled 
(massive) scalar. If the field is displaced from the minimum of its potential during inflation, 
there will be a phase where the field relaxes. Provided the scalar mass is much smaller than 
the curvature scale during inflation such a phase could be rather long. During the de Sitter 
stage, the specific form of the expanding background will dictate, through the equations of 
motion, the rate of suppression of the amplitude of 4>. 

While the field relaxes towards the minimum of its potential, energy is pumped from the 
homogeneous mode of (p to the gauge field fluctuations. The function /(</>) can be either 
an increasing function of (leading to a decreasing coupling) or a decreasing function of <f> 



1 The Heaviside electromagnetic system of units will be used throughout the investigation. The 
effective "electron" charge will then be given, in the present context, by e(<p) = ei/(0) -1 / 2 . 
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(leading to an increasing coupling). In both cases, depending upon the parameters of the 
model, the two-point correlation function of magnetic inhomogeneities increases during the 
inflationary stage. This implies that large scale magnetic fields can be potentially generated. 



Since the pioneering work of Fermi |15| large-scale magnetic fields are a crucial component 
of the interstellar and perhaps intergalactic medium |T7[. Faraday rotation measurements, 
Zeeman splitting estimates (when available) and synchrotron emission patterns conspire 
towards the conclusion that distant galaxies are endowed with a magnetic field of roughly 
the same strength of the one of the Milky Way 0,^8]. 

Observations of magnetic fields at even larger scales (i.e cluster, inter-cluster) are only 
at the beginning but, recently, promising measurements have been reported. Abell clusters 
with strong x-ray emission were studied using a twofold technique [19|,|2(|. From the ROSAT 
full sky survey the electron density along the line of sight has been determined. Faraday 
rotation (for the same set of 16 Abell clusters) has been determined through observations 
at the VLA []. The amusing result (confirming previous claims based only on one cluster 
|21|| ) is that x-ray bright Abell clusters possess a magnetic field of \i Gauss strength. The 
existence of strong magnetic fields with coherence scale larger than the galactic one can be 
of crucial importance for the propagation of high-energy cosmic rays. 

The large scale galactic magnetic fields are assumed to be the result of the exponential 



amplification (due to galactic rotation) of some primeval seed fields [T§j2^] . It was Harrison 
|T6| who suggested that these seeds might have something to do with cosmology in the same 
way as he suggested that the primordial spectrum of gravitational potential fluctuations 
(i.e. the Harrison-Zeldovich spectrum) might be produced in some primordial phase of the 



2 The ROetgen SATellite was flying from June 1991 to February 1999. ROSAT provided a map of 
the x-ray sky in the range 0.1-2.5 keV. 

3 The Very Large Array telescope is a radio-astronomical facility consisting of 27 parabolic antennas 
spread around 20 km in the New Mexico desert. 
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evolution of the Universe. Since then, several mechanisms have been invoked in order to 



explain the origin of these seeds |23],[^J and few of them are compatible with inflationary 
evolution. 

The plan of the present paper is then the following. In Section II the basic ideas con- 
cerning the model of evolution of the gauge coupling will be introduced. In Section III 
bounds coming both from the homogeneous and from the inhomogeneous evolution of will 
be described. In Section IV the evolution of the magnetic inhomogeneities will be addressed 
along the various stages of the model with particular attention to the role of the two-point 
function. The magnetohydrodynamical (MHD) approach will be generalized. In Section 
V the large scale magnetic fields produced in the scenario will be estimated. Section VI 
contains some concluding remarks. 



II. BASIC EQUATIONS 
A. Preliminaries 

Thanks to the high degree of isotropy and homogeneity of the observed Universe, the 
background geometry can be described using a (conformally flat) Friedmann-Robertson- 
Walker (FRW) line element 

ds 2 = G^dx^dx" = a 2 (ri)[dr] 2 - dx 2 }, (2.1) 

where r\ is the conformal time coordinate and G^ v is the four-dimensional space-time metric. 
The cosmic time coordinate (often employed in this investigation) is related to r] as a(rj) di] = 
dt. 

From the anisotropies of the Cosmic Microwave Background (CMB) it is consistent to 
assume that the Universe underwent a period of inflationary expansion of de Sitter or quasi- 
de Sitter type. Therefore a(r]) ~ —Vi/v during a phase stopping, approximately, when 
the curvature scale was Hi < lCC 6 Mp. For f] > f]i (possibly after a transient period) the 
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Universe gets dominated by radiation [i.e. a{rj) ~ 77] and then, after decoupling, by dust 
matter [i.e. a(rj) ~ rj 2 ]. 

The action describing the dynamics of the (Abelian) gauge coupling in a given back- 
ground geometry can be parametrised as 



S = d 4 x \[^G 



-G^d^d u 4> - V(cf>) - 



1 



-F UU F^ + J a A c 



(2.2) 



L2~ ~ pr ~" r ' yr/ 

In Eq. ( |2.2j ) the contribution of a classical (Ohmic) current has been included because of 
the relevance, for the problem at hand, of dissipative effects associated with the finite value 
of the conductivity, as it will be discussed in Section IV. 
From Eq. (|2.2|) the equations of motion can be derived 



1 



+ 



dV 



d<p 2g 3 ((j)) d(j) 
J?. 



dg F aB F«e + ^A« 



(2.3) 
(2.4) 



V^G a L<f(0) 

Since the gauge coupling appears in the Ohmic current, the term d$J a 7^ 0. Inserting Eq. 
( p.4| ) into Eq. ( |2.3| ) the explicit dependence upon the Ohmic current can be expressed, in 
terms of the gauge fields, as suggested in WM. Namely, using Eq. (|2.4p we can write that 



v^G^ = 2d u 



V^Gdg 



F m A r 



V^Gdg 



F a pF a ?. 



(2.5) 



An alternative form of Eq. ( |2.4|) becomes then 



7^ 



V^GG^d„ 



+ 



dV 



1 ^(0) vap 
2^(0) d(j) ■ 



(2.6) 



If the gauge coupling does not change, the evolution of Abelian gauge fields is conformally 
invariant. Hence, using the the conformal time coordinate, the appropriately rescaled gauge 
field amplitudes obey a set of equations which is exactly the one they would obey Minkowski 
space. In order to simplify the explicit form of the equations of motion the electric and 
magnetic fields will be rescaled in such a way that the obtained system of equations will 
reproduce the usual (conformally invariant) system in the limit g — > constant. The rescalings 
in the fields are 
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B = a 2 B, E = a 2 £, 



A = aA J = a 3 j, a = a c a, 



(2.7) 



where £>, £, A, j, a c are the flat-space quantities whereas B, E, A, J, a are the curved-space 



ones. 



Bearing in mind Eq. ( [2.7|) , the explicit form of Eqs. ( |3.10| )- (|2.6| ) becomes: 

1 dg 



d 2 cj) n M dV 2 - 
-4 + 2Hjf + —a 2 - V\ 
or] z or] ocp 



2g 3 a 2 



E 2 - B 2 



dB 
dr) 
d 



-V x E, 



dr] I g 2 (4>) 
V -5 = 



E 



J 



9 2 (<P) 



V x B--^-V<p x B 

gdcj) 



V-E = ~^E-V(f), 
gd(j) 

V-J = 0, J = a(E + vxB). 



(2.8) 

(2.9) 

(2.10) 
(2.11) 
(2.12) 
(2.13) 



where v is the bulk velocity of the plasma. The spatial gradients used in Eqs. ([2.8|)-( j2"TT3|) 
are defined according to the metric (|2.1|) . In Eq. ( p.!4|) the quantity 7i = dlna/dr] has also 
been introduced. 7i is the Hubble factor in conformal time which is related to the Hubble 
factor in cosmic time as7i = aH where H = dlna/dt. 

Once the background geometry is specified we are interested in the situation when the 
gauge field background is vanishing and the only fluctuations are the ones associate with 
the vacuum state of the Abelian gauge fields. Hence, Eqs. ( P . 8| ) ( p7T3[ ) allow to compute 
the evolution of cf> and the associated evolution of the two-point function of the gauge field 
fluctuations. 

Suppose that is originally displaced from the minimum of its potential. As far as the 
zero mode of <fi is concerned the system of equations can be further simplified: 



d 2 4> yyp4> 

drf dr] 



dV 



dg 



2g 3 a 2 ■ 



E — B 1 



dB 
dr] 



-V x E, 



(2.14) 
(2.15) 



d 



1 -E 



+ J = 4tvVx5, (2.16) 



drj ig 2 ((f)) J g 2 ((j)) 
V-5 = 0, V-£ = 0, (2.17) 

V-J = 0, J = a(E + vx B). (2.18) 

By now combining together the modified Maxwell's equations we obtain the evolution of the 
magnetic fields 

B" - 2—B' - V 2 -B = g 2 V x J, (2.19) 

where the prime denoted derivation with respect to the conformal time coordinate (the 
over-dot will denote, instead, derivation with respect to cosmic time). 

Once the evolution of the metric is specified, Eq. (|2.14j) dictates a specific evolution 



for and the evolution of <p will determine, in its turn, the evolution of the gauge fields. 
The interesting initial conditions for the system are the ones where the classical gauge field 
background vanishes. Thus, when the homogeneous component of <p starts its evolution 
during the de Sitter phase, quantum mechanical fluctuations will be postulated as initial 
conditions of gauge inhomogeneities. 

When the background geometry evolves from the de Sitter phase to the subsequent epoch, 
massive quanta of are produced. The amount of the produced inhomogeneous modes of 
cf) can be computed and it will be shown that the associated energy density will always be 
smaller than the one of the homogeneous mode. This analysis will be one of the subjects 
discussed in Section III. 



III. CONSTRAINTS ON THE EVOLUTION OF THE GAUGE COUPLING 

Sub-millimiter tests of the Newton's law show that no deviations are observed down to 
distances as small as 0.1 mm [p6|| . Therefore, <fi should be massive. Of course the potential 
of (ft may be much more complicated than the one provided by a simple mass term. However, 
for sake of simplicity, a massive scalar will be analyzed since, already in this case, interesting 
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effects can be analyzed. In spite of the fact that this choice is apparently simple, various 
constraints on the scalar mass appear. 



A. Evolution of the homogeneous mode 

Suppose that the potential term driving the evolution of the gauge coupling is simply 

2 

77? 

V{<t>) ~ — 2 . (3.1) 

During the inflationary stage of expansion the scale factor evolves as a(rj) = (—771/77) for 
77 < —r] 1 . The evolution of (ft is obtained by solving 

2 a 2 

0" + -0' + ^0 = 0, (3.2) 

7/ T] z 

where \x = m/H\ and where the relation Ti ~ r/ _1 ~ aH has been used. If // <C 1, for 



77 < —7/1 the solution of Eq. Q3.2|) can be written as 



i (7/)=0 1 -0 2 (-^-) . (3.3) 

The end of the inflationary stage of expansion may not be directly followed by the radiation 
dominated phase. In the intermediate phase the scalar mass is still small than the curvature 
scale but the curvature decreases, in general, faster than during the inflationary phase since 
the background is neither of de Sitter nor of quasi-de Sitter type. The evolution of the scale 
factor can be parametrised as 0(77) ~ 77° where, in order to fix the ideas, a ~ 2 could be 
assumed 0. The evolution of will simply be 

\la — 4 3 /7?\ 2a_1 l 

4> A (V) = 01 + 02 T + ~ " - , T]l<V<Vr, (3.4) 

Via — 1 la — 1 \rjiJ J 
where the continuity between Eq. p.3j ) and Eq. Q3.4Q has been required, so that 0i(— r/x) = 
rh (7?i) and (f)[(-r)i) = <j)' lh (vi)- 



l The case a = 2 corresponds to a matter-dominated intermediate stage. 
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After 7] r the background enters a radiation dominated phase and the evolution of can 
be explicitly solved in cosmic time. The equation for <p , in this phase, is given by 



+ 3#0 + m 2 = 0, H 



(3.5) 



which in terms $ = a 2 0, becomes 



m 2 _ _ 
2 4 



$ = 0. 



(3.6) 



In the radiation-dominated stage of expansion Eq. ( |3.6| ) becomes 

3 



m 2 + 



16t 2 



$ = 0, 



(3.7) 



whose solution can be written in terms of Bessel functions [ 27fl 



$(m£) 



AYi (mt) + 5 Ji (mt) 



(3.8) 



For mt <C 1, 4> has a constant mode and a solution as t~ 1//2 . Recalling the relation between 
cosmic and conformal time and imposing the continuity of and 0' (in rj r ) with the solution 
of Eq. (|3.4| ) the following form can be obtained: 



+ 



2a -4 6(1 - a) fr]^ 2 ^ 1 



(3.9) 



A 2a — 1 2a — 1 \r] r J J Vr/ r / 77 J 

which is valid for r] r < 77 < r/ m . The time ^ m marks the moment where if ~ m. When 
mt > 1, the regime of coherent oscillations takes over and the solution (|3.8|) implies that the 
energy density stored in decreases as a -3 , meaning that (p c (v) ~ i]~ 3 ^ 2 . Since the coherent 
oscillations decrease as a -3 there will be a typical curvature scale H c and a typical time r] c at 
which the coherent oscillations become dominant with respect to the radiation background. 
This moment is determined by demanding that 



^M^(^ 4 ~m 2 2 f^ 3 
\a c / V a c 



(3.10) 



which also implies that 



H c ~ rrup 
10 



(3.11) 



where ip = 0i/Mp. Eq. (|3.11|) has been obtained without tuning the asymptotic value of 
cf) to the minimum of its potential. If such a tuning is made, the amplitude of oscillations 
at r/ m will be smaller than <pi and it will be given, according to Eq. ( |3.9| ), by 4>2(vi/ r tr) 2a l ■ 
Thus, the scale H c will be defined by a different relation namely: 

leading, ultimately, to 

L A TT 4(2q-l) 

) 2 \ 4 / H r \ a+i / m 



In the approximation of instantaneous reheating [i.e. 1] V ~ rji], H r ~ H\. Therefore, from 
Eq. (|3.13|) H c is smaller than the value determined in Eq. ( |3.11| ) by a factor (m/Hi) 2 . In 



the approximation of matter-dominated reheating (i.e. a ~ 2), the result of instantaneous 
reheating is further suppressed by a factor (H T / Hi) 4 as one can easily argue from Eq. Q3.13|) . 



From Eqs. ( |3.5| )-( |3l)D , the evolution of <fi will go as r/ -3 when coherent oscillations start 
dominating. 

In spite of the possible tunings made in the asymptotic values of 0, after rj c there will 
be a typical time at which the field <fi will decay. In order not to spoil the light elements 
produced at the epoch of BBN has to decay at a scale larger than if ns ~ T^/Mp (where 
T ns ~ MeV. Since <fi is only coupled gravitationally the typical decay scale will be given by 
comparing the rate with the curvature scale giving that 

in 

Ht~r~ — >Hns , (3.14) 

implying that m > 10 4 GeV. This requirement also demands that the reheating temperature 
associated with the decay of <fi will be larger than the BBN temperature. 

In order to illustrate some concrete examples of the various possibilities implied by our 
considerations, suppose that m ~ 10 3 TeV and suppose that the asymptotic value of <ft is 
fine-tuned to its minimum. Furthermore, suppose that the reheating is instantaneous. Then, 
according to the picture which has been presented, inflation stops at a scale Hi ~ 10 13 GeV 
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and starts oscillating at a curvature scale H m ~ 10 3 TeV. The coherent oscillations will 
then become dominant at a curvature scale H c ~ 10~ 8 GeV (having assumed 02 ~ A4p). The 
coherent oscillations of will last down to H$ ~ 10~ 20 GeV. After this moment the Universe 
will be dominated by the radiation produced in the decay of 0. Notice, for comparison, that 
the BBN curvature scale is H ns ~ 10~ 25 GeV so that the decay occurs well before BBN (five 
orders of magnitude in curvature scale). 

Another illustrative example is the one where m ~ 10 6 TeV. In this case the decay of 
occurs prior to the EWPT epoch, namely 

H 4> >H ew . (3.15) 

In fact H ew = v^V^T 2 w /M P ~ 1(T 17 GeV (with N eS = 106.75 and T ew ~ 100 GeV) whereas, 
from Eq. ( |3.14j) , ~ 10~ 9 GeV. In more general terms we can say that in order to have 



the decay occurring prior to the EWPT epoch we have to demand that > H ew which 
means that m > 10 5 TeV. 

In closing this section two general comments are in order. If no fine-tuning is made in the 
asymptotic amplitude of 0, the typical scale of the coherent oscillations will almost coincide 
with m. However, the possibility tp <C 1 is still left if, for some reason, we want H c m. 

The decay of and the consequent freezing of the gauge coupling should occur prior 
to the EWPT epoch and the baryon number should be generated, in the present context, 
at the electroweak time [§,0. Suppose, for example, that this is not the case and that the 
BAU has been created prior to the electroweak scale. Suppose, moreover, that the decay of 
occurs after baryogenesis. Then the temperature of the radiation gas before the decay of 
4> will be ~ T m (a m /a^) ~ m(m/Mp) 5 . Thus, the entropy increase due to the decay of 



will be AS ~ (^decay/^) 3 where T decay ~ ^jH^M-p. This implies that AS ~ m/Mp. It has 
been observed in different contexts that in order to preserve a pre-existing BAU one should 
have AS < 10 5 |29| , |30| . Thus, this bound would imply m > 10 14 GeV. This is the reason 
why the present analysis will assume that the decay of occurs prior to the electroweak 
time and that the BAU is generated at the EWPT or shortly after. 
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B. Evolution of the inhomogeneous modes 



When the Universe passes from the inflationary stage to the subsequent radiation dom- 
inated expansion, inhomogeneities of the field <fi are generated. This may invalidate the 
original assumptions and introduce further complications by adding qualitatively new con- 
straints on the scenario. 

It is useful to recall that the inhomogeneities of 4> can be interpreted, in the framework of 
second quantization, as quanta of the field 0. Hence, the inhomogeneities produced because 
of the sudden change of the geometry from the de Sitter epoch to the radiation dominated 
epoch can be counted by estimating the number of quanta produced by the sudden change 



of the geometry according to the well known techniques of curved space-times [[31 
Consider the first order fluctuations of the field 

<f>{x,ri) = <j>[rj) + 5(j)(x,T]), (3.16) 

whose evolution equation is, in Fourier space, 

i)" + 2Hip' + [k 2 + m V]^ = 0, (3.17) 

where ip(k, rj) is the Fourier component of 5<p(x, rj). In order to count the number of quanta 
produced during the transition of the geometry from the inflationary to the radiation dom- 
inated stage of expansion the (canonically normalized) amplitude of fluctuations \I/ = ipa 
should be defined so that Eq. (|3.17f) becomes: 



a" 



f +[fc 2 + mY--]$ = 0. (3.18) 
a 



In the de Sitter stage of expansion Eq. fl3.18|) reduces to 

^ + [k 2 + ^-^)% = 0, (3.19) 
whereas during the radiation dominated stage of expansion Eq. (|3.19| ) takes the form 

Vi 
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The solution of Eq. ( |3.19| ) (with the correct quantum-mechanical normalization for rj — > 
— oo) can be written as 

tfifo) = -L V v 7 ^ Hf\-x), (3.21) 
where x = krj and is the first order Hankel function p7fl . In the pure de Sitter case, 



p = 3/2-y/l — (4/9)/i 2 and since /x<Cl,p — 3/2;pisa phase factor which has been chosen 
in such a way that 

p= y| e *5d+2p). (3.22) 

With this choice of p we have that ^(77) ~ e~ lkv / y/2k for 77 — > —00. 

During the radiation dominated stage of expansion Eq. (3.20) is the equation of parabolic 



cylinder functions [27]. The solutions turning into positive and negative frequencies for 
77 — > +00 are then 

f*(v) = (^171 e " f D ig ^(e^z), (3.23) 

where 

z=^(r ] + 2r ]l ), q=^ (3.24) 

27 

and where D a are the parabolic cylinder functions in the Whittaker's notation. The solution 
ofEq. (|3~20| ) 

^(77) = c + (% r (r7) + a-(k)g*(ri), (3.25) 

is given in terms of c+(k) and C-(k) which are the two (complex) Bogoliubov coefficients 
satisfying |c + (/c)| 2 — |c_| 2 = 1. In a second quantized approach |c_(fc)| 2 is the mean number 
of created quanta, whereas in a semi-classical approach c_(/c) can be viewed as the coeffi- 
cient parametrising the mixing between positive and negative frequency modes. In the case 
c_(&) ~ no mixing takes place and no amplification is produced. In order to determine 
c±(k), ^i(rj) and \I/ r (77) should be continuously matched in 77 = —771, namely 
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(3.26) 



By solving this system, an exact expression for the Bogoliubov coefficients is obtained which 
is, in general a function of two variables : p = mr]i and x% = krji. Since /i < 1 the exact 
result can be expanded, in this limit, 

c+ (k) = 7T e i %{-^J^S 2 (x u p)fi- 1 * + ( r^jy[Si{x 1 , p) + S 2 (x!, p)]fi*} + 0(ji*), 

c _(k) = ne^{--jJ^S,(x u p)vrt + ^W 1)P ) + S 2 (x 1 , p)]p} + O(pi), (3.27) 

where Si(xi, p) and 5*2 (xi, p) contain the explicit dependence upon the Hankel's functions: 



S 1 (x 1 ,p) = e i ^ 1+ ^H^(x 1 ), 



S 2 (x 1 ,p) = y ^ 1 e i ^ + ^\(p+l) H ^ {Xl] 



x 1 H i p_l 1 (x 1 



(3.28) 



If p ~ 3/2 Eq. (]pD gives 



c+(k) = e*8 Vtt 
c-(k) = e~ l ^\pK 



x 



- + 



% X 



- + 



1 (i - 1W/ 4 

r + 



2Y@pi 2r(f)/i4 L 2 r(|)^ V2T{\ 

■ 1 % X i 

T + 



(3.29) 



2r(|)^ 2r(|)^ L 2r(J)A*4 v^r® 

In the limit x% = krji <C 1 the mean number of created quanta can be finally approximated 



as 



n(k) ~ |c-(fc)| =g|A;r/i| p /i 



-2p -1/2 



(3.30) 



where g is a numerical coefficient of the order of 10 2 . The energy density of the created 
(massive) quanta can be estimated from 

d 3 u 



dp^ 



(2vr) : 



m n(k) 



(3.31) 



where u = k/a is the physical momentum. In the case of a de Sitter phase (p = 3/2) the 
typical energy density of the produced fluctuations is 
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(^gmfff^)^^) 3 (3.32) 



The produced massive quanta may become dominant. If they become dominant after 
already decayed they will not lead to further constraints on the scenario. If they become 
dominant prior to the decay of <fi further constraints may be envisaged. The scale at which 
the massive fluctuations become dominant with respect to the radiation background can be 
determined by requiring that p^(?7*) — Piiv*) implying that 

which translates into 

H*~q 2 me 4 , (3.34) 

where e = H\/Mp. In order to make sure that the non-relativistic modes will become 
dominant after already decayed < should be imposed, that is to say m > 10 2 TeV 
for e ~ 1(T 6 . 

The maximum tolerable amount of entropy, in order not to wash-out any preexisting BAU 
is model-dependent but, in general, AS < 10 5 seems to be acceptable p8|-p0[]. Defining 

as the radiation gas already present at the scale H^, the entropy increase from to 
T decay ~ ^H^M-p is of the order of 



T \ 

- 1 decay \ 



AS = > (3-35) 
where 

T, = tJ^) ~ m e /6 (3.36) 
where £ = m/M-p. Demanding that AS < 10 5 implies that 

£ > i(T 10 e 3 . (3.37) 
Taking, as usual, e ~ 1CT 6 , m > 10 GeV. 
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Hence, if the constraints pertaining to the homogeneous mode are enforced, the bounds 
coming from the inhomogeneous modes do not invalidate the conclusions of the analysis. 



According to the logic expressed in Eq. (|3.15 ), an illustrative example is the case where 
m > 10 5 TeV and the BAU is generated after EWPT. In this case the bounds obtained in the 
present section are satisfied and the analysis of the evolution of the inhomogeneous modes 
shows that the qualitatively new bounds introduced in the picture are less constraining than 
the ones obtained in the analysis of the dynamics of the homogeneous mode. 



IV. EVOLUTION OF THE GAUGE FIELD FLUCTUATIONS 



The evolution of the field <ft during and after the de Sitter stage implies, according to 



Eqs. ( |3.3| )-( f01 ), that the two-point function of the gauge field fluctuations may very well 
grow. 

For r) < —r)i, a rough approximation suggests that the effect of the ohmic current is not 
present and the gauge field is in the vacuum state with k/2 energy in each of its modes. By 
promoting the classical fields to quantum mechanical operators we have that the physical 
polarizations of the magnetic field can be written as 

/H 3 k 
<[a k , a b(k V )e^ + a\j?(krj)e- ik -% (4.1) 



where b(krj) = B{krj) / g(rj) obey the equation 



b" + 



9 J 9 J 



6 = 0. (4.2) 



Notice that b(krj) are the correct normal modes whose limit (for rj — > — oo) should be 



normalized to ^k/2e~ lkv . The two-point correlation function of the magnetic fluctuations 
can then be expressed as 

r AZh 

Q^n) = Ux^^x + r^)) = j Pij (k)b(k, V )b*(k, V )e ik \ (4.3) 



where 
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PM = (Si, - (4.4) 

The magnetic energy density, derived from the energy-momentum tensor corresponding to 
the action of Eq. Q2.2p , is related to the trace of the correlation function reported in Eq. 
Qi.lD over the physical polarizations: 

/sm J\iir dk 
PB(k,T})— — - (4.5) 

where 

p B (k, V ) = ^k 3 \b(k,7 1 )\ 2 . (4.6) 

7T 

A necessary condition in order to assess that gauge field fluctuations grow during the de 
Sitter phase is that the two-point function increases in the limit i] — > — rji 



Suppose that the gauge coupling decreases with monotonic dependence upon the field </>, 
namely 

This parametrisation is purely phenomenological, however, it allows to take into account, 
at once, some physically interesting cases like the one suggested by the low-energy string 
effective action where, in the limit of 0/Mp < 1, g 2 (4>) ~ <p. 

Using Eq. fl3.3| ) and Eq. ( |4.7| ) into Eq. ( |4.2| ) the time evolution of the normal modes of 
the magnetic field can be found analytically since the specific form of Eq. ( |4.2| ) falls in the 
same category of Eq. ( |3 . 1 9| ) . Hence, 



b(k,ri) = Ny/krjHjp(kq), N = X^ e -^ (1+2i/) , (4.8) 



with v = (3A + 1J/2 and where Hi (kr)) the Hankel function of second kind |p7|1 . Notice that 
the normalization N has been chosen in such a way that for rj —>■ — oo the correct quantum 
mechanical normalization is reproduced. Consequently, following Eq. ( |4.3| ), the two-point 
function evolves as 
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— 3A 

lim 0tf(r,tj) ~ . (4.9) 

Since the two-point function increases magnetic fluctuations are generated. 

For sake of completeness the case of increasing gauge coupling will now be examined 
using the following phenomenological parameterisation 

Again different scenarios can be imagined. For instance, one coould argue in favour of 
scenarios where the gauge coupling depends upon <fi (or upon rf) in a highly non-monotonic 
way. For the illustrative purposes of the present investigation it is however sufficient to focus 
the attention on the case of monotonic dependence. 

Following now the same steps outlined in the case of decreasing gauge coupling the 
evolution of the two-point function can be obtained 

rt-i 2 — 36 

lim Gij{r,rj) ~ ± . (4.11) 



n-*-vi 



for S > 1/3 and 



77 3S 

lim Gij(T,-n)~ . (4.12) 

n-t-rn rji 

for 5 < 1/3. If 8 < 1/3 the correlation function decreases and this signals that large scale 
magnetic fields are not produced. 

The back-reaction of the produced fluctuations can be safely neglected in de Sitter space. 
Looking at Eqs. @ and (|2lj) it can happen that if the magnetic fluctuations grow too 



much the term at the right hand side will become of the same order of the others. This is 
not the case. Using the conventions of this Section together with the explicit form of the 
scale factor in the de Sitter phase it can be shown that 

-2v 



1 d9B2 ' \k Vl \ 5 ~ 2l/ (4.13) 



g 3 a 2 



-B- ~ 



where v is determined from the specific power dependence of the coupling as a function of <fi. 
Since we are interested in large scale modes we have kr]i 1. Therefore the back reaction 
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effects are relevant towards the end of the de Sitter phase (i.e. r] ~ — r/i ) and for k ~ r)i , 
namely exactly for the modes not relevant for the present investigation. 

After the end of inflation the onset of the conductivity dominated regime may not be 
instantaneous. In this case after rji the presence of a reheating phase should be taken 



into account. Suppose, for instance, that g{rj) decreases according to Eq. (|4.7|). Suppose, 
moreover, that during reheating the background is dominated by the coherent oscillations 
of the inflaton. In this case the effective evolution of the geometry will be dominated by 
matter with a(r]) ~ r] 2 . According to Eq. fl3.4|) (with a ~ 2), ~ r]~ 3 . The value of the 
magnetic inhomogeneities at r/ r will be given by solving Eq. ( |4.2|) 



3 



A / n \ 1-#A 



b(k,r} T ) ~ AA^Y +A 2 (^) - , (4.14) 



,771/ 

where Ai and A 2 are two arbitrary constants. Depending upon the value of A the fastest 
growing solution is selected in Eq. (|4.14j ). This phase may induce further amplification on 
the two-point function since its main effect is to delay the conductivity-dominated regime. 



A. Generalized MHD equations 

For 1] > r] r the role of the Ohmic diffusion becomes important and the evolution of the 
magnetic inhomogeneities will be described by the MHD equations generalized to the case 
of time varying gauge coupling. 

The ordinary (i.e. fixed coupling) MHD treatment is an effective description valid for 
length scales larger than the Debye radius and for frequencies smaller than the iono-acoustic 
frequency. This means that MHD is accurate in reproducing the spectrum of plasma excita- 
tions obtained from the full kinetic (Vlasov-Landau) approach but only for sufficiently low 
frequencies and for sufficiently large scales. 

Implicit in the ordinary MHD analysis is the assumption that the plasma has to be 
electrically neutral (V • E — 0) over length scales larger than the Debye radius. Thus, 
this system of equations cannot be applied for distances shorter than the Debye radius and 
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for frequencies larger than the plasma frequency where a kinetic description should be 
employed. 

MHD equations can be derived from a microscopic (kinetic) approach and also from a 
macroscopic approach where the displacement current is neglected ||33|| . If the displacement 
current is neglected the electric field can be expressed using the Ohm law and the magnetic 
diffusivity equation can be derived 

^ = Vx(?;x5) + -V 2 5. (4.15) 
or] a 

The term containing the bulk velocity field is called dynamo term and it receives contribution 
provided parity is globally broken over the physical size of the plasma. In Eq. (|4.15|) the 



contribution containing the conductivity is usually called magnetic diffusivity term. 

In the superconducting (or ideal) approximation the resistivity of the plasma goes to zero 
and the induced (Ohmic) electric field is orthogonal both to the bulk velocity of the plasma 
and to the magnetic field [i. e. E ~ — v x B]. In the real (or resistive) approximation the 
resistivity may be very small but it is always finite and the Ohmic field can be expressed as 

E ~ V X B -vx B. (4.16) 
a 

If the gauge coupling changes with time the system of equations obtained by neglect- 
ing the displacement current receives new contributions and the relevant equations can be 
obtained, in the resistive approximation, from Eqs. ( 2.14|) -( ^.18 ): 



VxE=-^- ) (4.17) 

E = --vxB, (4.18) 
a 

1 - q 
— V x B = 



2 VxB = J-2^E, (4.19) 



Using Eqs. ( |4.17| )-( [4.19j ) the generalized magnetic diffusivity equation can be obtained: 

(1 - -~)%r = Vx(vxB) + —,V 2 B. (4.20) 
a g 2 g dr] ag 2 

Notice that Eq. fl£20p reproduces Eq. fl£l5|) if g' -»• 0. 
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Suppose now that the plasma, whose effective Ohmic description has been presented, is 
relativistic. In the case when the coupling is constant a is constant and it is given by 

a = a c ( V )a( V ), (4.21) 

where a c ~ T/g 2 scales as the inverse of a(rj) if the evolution of the Universe is, to a good 
approximation, adiabatic. 

If g is not constant, a is not constant anymore but it decreases if the gauge coupling 
increases and, vice versa, it increases if the gauge coupling decreases. In spite of this, in the 
generalized MHD equations, the combination which appears is always ag 2 which is roughly 
constant for an adiabatically expanding Universe. 

In the approximation of instantaneous reheating, the solution of Eq. ( |4.20| ) is given by 

-J fc2 , drj 

B(k,T)) = B{k,r]i)e ° g2 - 2a v . (4.22) 



According to Eqs. (|3.13| )- (|3.14j) , in order to get the coupling frozen prior to H ew ~ 1CT 17 
GeV, m > 10 5 TeV shall be required. If the gauge coupling is always decreasing as a function 
of rj, it can be parametrized by Eq. ( |4.7| ). Hence Eq. ( ^4.22| ) can be evaluated by using the 
explicit evolution of 4> as obtained from Eqs. ( |3 . S| ) — ( PTD|) implying that (fi T ~ 77 -1 , m ~ ry~ 3 / 2 
and C ~ r]~ 3 . The result is 

B(k, Vo ) =I{ Vl , 77 m ,7fo, Vc )e-^ {m+ ' no) (4.23) 

where 770 is the present time and 

A + ag 2 i] m 3A + 2ag 2 i] c 3A + ag 1 ^ 



Concerning Eqs. ( [4.23| )-( |4.24| ) few comments are in order. From Eq. ( [4.23| ) all the modes 



k 2 > k 2 a ~ ^ (4.25) 
Vo 

are suppressed by the effect of the conductivity. The present value of uv(?7o) can be 



5 With Ld(rj) ~ k/a(rj) the physical momentum will be denoted. 
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estimated by recalling that I/770 ~ Hqcio where Hq ~ 10~ 61 Mp. Thus uo a ~ 10~ 3 Hz. 
Present modes of the magnetic fields are dissipated if to > u a . 

As far as the problem of galactic magnetic fields is concerned, the relevant set of scales 
range around the Mpc corresponding to present modes of the magnetic field ujg ~ 10~ 14 Hz, 
i.e. uj g <ti uj a . 



V. MAGNETIC FIELD GENERATION 

Large scale magnetic fields have been postulated over fifty years ago in the context of 
cosmic rays. Prior to the development of magnetohydrodynamics there was the common 



belief that cosmic rays are in equilibrium with the stars (like the sun) |15]. Fermi was 
probably the first one to realize that if our galaxy has a magnetic field, then cosmic rays 
could be in equilibrium with the large scale magnetic field of the galaxy. Implicit in the 
Fermi argument there was also the idea that large scale magnetic fields could also be present 



in other galaxies and, indeed, Fermi and Chandrasekar tried to develop the first theory 
of gravitational instability in the presence of large scale magnetic fileds. Today large scale 
magnetic fields are measured with a number of different techniques. Among them Faraday 
rotation is certainly the most robust. 



A. Generalities 

When polarized radiation passes through a cold plasma containing a magnetic field the 
polarization plane of the incoming radiation gets rotated by an amount which is directly 
proportional to the Larmor frequency (and then to the magnetic field), to the square of the 
plasma frequency (and then to the electron density). The quantity which is experimentally 
measured is the shift in the polarization plane per quadratic interval of wave-length, namely: 

RM = %L = 811.9 / f-^-W T 4 (5.D 

AA^ J Vcm 6 J VyuGauss/ m 2 
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where £ = Z/kpc and I is the integration variable running over the line of sight. It is 
clear from the expression of Eq. ( |5.1| ) that an independent estimate of the electron density 
is needed in order to asses the precise value of the magnetic fields we ought to measure. 
In the interstellar medium an estimate of the electron density along the line of sight is 
provided by the so called dispersion measurements JT7| , [TE| ]. Pulsars emit regular pulses of 
electromagnetic radiation with periods ranging from few seconds to few milliseconds. By 
comparing the arrival times of different radio-signals at different radio-wavelengths, it is 
found that they are slightly delayed as they pass through the interstellar medium exactly 
because electromagnetic waves travel faster in the vacuum than in an ionized medium. From 
the amount of the delay one can infer the dispersion measurement , namely, j n e dl . Dividing 
the rotation measurement by the dispersion measurement the mean magnetic field along the 
line of sight can be estimated. 



Recently magnetic fields have been measured in the intra-cluster medium ||19|| . ROSAT 
satellite identified a number of x-rays bright Abell clusters (XRBAC). From the surface 
brightness of the cluster the thermal electron density can be obtained. Various XRBAC 
(sixteen) have been monitored also with VLA observations. The clusters have been selected 
in order to show similar morphological features. From VLA observations the RM has been 
obtained, and from ROSAT the electron density has been obtained. The results show the 
presence of large scale magnetic fields ( of /xGauss strength) in the intra-cluster medium. 
The possible existence of large scale magnetic fields beyond the galaxy is also rather crucial 
for the deflection of high energy cosmic rays |35|j3^ • 



B. Origin of large-scale magnetic fields 

The discovery of large scale magnetic fields in the intra-cluster medium implies some 
interesting problems for the mechanisms of generation of large scale magnetic fields. Let us 
consider, first of all, magnetic fields in galaxies. Usually the picture for the formation of 
galactic magnetic fields is related to the possibility of implementing the dynamo mechanism. 
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The galaxy rotates with a typical period of 10 8 yrs [ I7[ . By comparing the rotation period 
with the age of the galaxy (for a Universe with ~ 0.7, h ~ 0.65 and Q m ~ 0.3) the 
number of rotations performed by the galaxy since its origin is approximately 30. During 
these 30 rotations the dynamo term of Eq. ( |4 . 1 5| ) dominates against the magnetic diffusivity 
term since parity is globally broken over the physical size of the galaxy. As a consequence 
an instability develops. This instability can be used in order to drive the magnetic field 
from some small initial condition up to its observed value. Most of the work in the context 
of the dynamo theory focuses on reproducing the correct features of the magnetic field of 
our galaxy |L7|,[18|. For instance one could ask the dynamo codes to reproduce the specific 
ratio between the poloidal and toroidal amplitudes of the magnetic field of the Milky Way. 

In spite of these aspects, if the number of rotations of the galaxy is approximately 30, 
the achievable amplification produced by the dynamo instability can be at most of 10 13 , i.e. 
e 30 . Thus, if the present value of the galactic magnetic field is 10~ 6 Gauss, its value right 
after the gravitational collapse of the protogalaxy might have been as small as 10 -19 Gauss 
over a typical scale of 30-100 kpc. 

There is a simple way to relate the value of the magnetic fields right after gravitational 
collapse to the value of the magnetic field right before gravitational collapse. Since the 
gravitational collapse occurs at high conductivity the magnetic flux and the magnetic helicity 
are both conserved. Right before the formation of the galaxy a patch of matter of roughly 1 
Mpc collapses by gravitational instability. Right before the collapse the mean energy density 
of the patch, stored in matter, is of the order of the critical density of the Universe. Right 
after collapse the mean matter density of the protogalaxy is, approximately, six orders of 
magnitude larger than the critical density. 

Since the physical size of the patch decreases from 1 Mpc to 30 kpc the magnetic field 
increases, because of flux conservation, of a factor (p a /Pb) 2 ^ 3 ~ 10 4 where p a and pb are, 
respectively the energy densities right after and right before gravitational collapse. Hence- 
forth, the correct initial condition in order to turn on the dynamo instability is B ~ 10 -23 
Gauss over a scale of 1 Mpc, right before gravitational collapse. 
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Since the flux is conserved the ratio between the physical magnetic energy density and 
the energy density sitting in radiation is almost constant and therefore, in terms of this 
quantity (which is only scale dependent but not time dependent), the dynamo requirement 
can be rephrased as 

tb{L) = E^L^p. > io~ 34 , L~lMpc. (5.2) 

P-rW 

If the dynamo is not invoked but the galactic magnetic field directly generated through some 
mechanism the correct value to impose at the onset of gravitational collapse is 

r B (L)~10" 8 , L~lMpc. (5.3) 



Clearly even the number given in Eq. Q5.2|) is, physically, not so small. The magnetic energy 
density stored in a quantum mechanical fluctuations of 1 Mpc is, in terms of re, 1CT 95 . 

The possible applications of dynamo mechanism to clusters is still under debate and it 
seems more problematic [[37]]. The typical scale of the gravitational collapse of a cluster is 



larger (roughly by one order of magnitude) than the scale of gravitational collapse of the 
protogalaxy. Furthermore, the mean mass density within the Abell radius ( ~ l.hh~ l Mpc) 
is roughly 10 3 larger than the critical density |55fl . Consequently, clusters rotate less than 
galaxies since their origin and the value of t-b(L) has to be larger than in the case of galaxies. 
Since the details of the dynamo mechanism applied to clusters are not clear, at present, it 
will be required that t-q{L) ^> 10~ 34 (for instance r^(L) ~ 10~ 20 ), in order to see if in the 
parameter space of the present model magnetic fields larger than the (galactic) dynamo 
requirement. 

C. Estimates of large scale magnetic fields 

In the approximation of instantaneous reheating the typical (present) frequency corre- 
sponding to the end of inflation can be computed and it turns out to be 

^oH^T^e 1 / 2 ^-!. (5 .4) 
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Since Td ec ~ 0.26 eV , z^T^ec ~ 100 GHz. Eq. (|5.4j ) can be obtained by red-shifting the 
highest mode, i.e. ^i(?7i) ~ Hi through the different stages of the evolution of the model, 
namely, according to Eqs. (|3.11| ) and (|3.14|) , from T]i down to r] m and from rj m down to r] rmc . 
Recall that from r] c to rj^ the Universe is, effectively, matter dominated. The other typical 
frequencies appearing in the time evolution of the gauge coupling can be written, in units 
of ux(r)o), as 

u m(Vo) _ _i/2 £-1/2 



Uc(Vo) {1/2 -1/2 



] £ -l/2 £ 7 /6 ^ (55) 



u <j>(V0) _ -1/2 t7/6 2/3 



Wi(?/o) 

where, as in the case of Eq. ( |5.4| ) all the frequencies are evaluated at the present time. 
In the case of decreasing gauge coupling [described by Eq. ([4.71)1 the amount of generated 



large scale magnetic field can be estimated from Eqs. ( |4.2| )-( [18D together with Eqs. ( [4.20| )- 
(fOS). Bearing in mind that the typical frequency scale corresponding to 1 Mpc is 10~ 14 



Hz we have, following the notations of Eq. (|5.2|) 

t b (uj g ) = /(A) 10-^-^ T{u G ), (5.6) 



where 



and 



Q3A-1 o 1 \ 



r(wG ) = e4[(^)(«*)i(^)§r4, (5.8) 



which means, using Eqs. (|5.4|)-(p3|). 



T(u G ) = e~^[e- l ' 2 V - i e 12 } (5.9) 



Notice that in Eq. (5J3) is the present CMB temperature. 
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In order to illustrate the regions of the parameter space where magnetogenesis is possible 
<f ~ 1 will be assumed. In Fig. [I] and Fig. ^ the case of decreasing gauge coupling is 
discussed. In Fig. [I], A is fixed and the exclusion plot is given in terms of e and £. In 
particular, for illustration, A = 1 has been chosen. The choice of A ~ 1 implies that 
g 2 (4>) ~ 0. Such a case is favoured from the tree- level string effective action where the 
effective coupling can be approximated as g 2 (4>) ~ 0/Mp for < Mp. 

The two vertical lines lines mark, respectively, the bounds coming from BBN [i.e. £ > 
and from the electroweak epoch [i.e. £ > 10 , from Eq. (|3.15| )1. 
For consistency with the assumptions of Eqs. (|3.2j )- (|3.3| ) m/Hi <C 1, implying e ^> £. With 
the lower dashed line the bound e > £ is reported. Notice that the requirement of Eq. (|3.37|) 
is not numerically relevant. The upper dashed line is obtained by requiring, according to 
Eq. ( |3.34| ), that > H*. Finally, the full (diagonal) line is derived by imposing on Eqs. 
fl5.6|)-(jOD the dynamo requirement of Eq. (|5.2|). 

In Fig. [I the value of £ has been fixed to 10~ n , as required by the considerations related 
to the electroweak epoch and the magnetogenesis region is described in terms of e and A, 
both varying over their physical range. The two horizontal lines fix the bounds coming from 
e < 10~ 6 and from e > £. With the dashed line the curve r s (c<j G ) ~ 10~ 20 is denoted. Hence, 
values larger than the dynamo requirement are allowed. This observation may be relevant 
in the context of magnetic fields associated with clusters. In Fig. ^|, A lies in the range 
< A < 4/3. This choice guarantees the growth of the correlation function of the magnetic 
inhomogeneities during the de Sitter stage. 

In order not to conflict with large scale bounds coming from the isotropy of the CMB the 
energy spectra of the produced gauge field fluctuations have to decay at large distance scales, 
implying that < A < 4/3. To be compatible with CMB anisotropies r^(u>dec) < 10~ 10 
should be imposed (where c^dcc — 10 -16 Hz). If the condition < A < 4/3 is enforced, the 
spectra increase with frequency and the possible bounds coming from the anisotropy of the 
CMB are satisfied. 

An analogous estimate can be obtained in the case of increasing gauge coupling discussed 
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10~ 15 , from Eq. (3.14 
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FIG. 1. The shaded area illustrates the region where magnetogenesis is possible in the case 
where A = 1 and (/? ~ 1. The vertical lines correspond to the requirements coming from BBN and 
from the EWPT epoch. 
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FIG. 2. The region with stripes defines the area where magnetogenesis can occur in the case of 
fixed mass [i.e. £ ~ 10~ n ] and for tp ~ 1. 



in Eq. Q4.10 ). In this case 



(5.10) 



where 



f(S) 



7T 3 V2 2 



(5.11) 



and 



2 



(5.12) 



In Fig. |3| and Fig. [| the requirements coming from the dynamo mechanism as well as the 
other theoretical constraints are illustrated. For both plots e < 10~ 6 and ip ~ 1. In Fig. 
H] the case 5 = 2 is illustrated. The shaded area selects the region of the parameter space 
where the dynamo requirement of Eq. (|5.2| ) is imposed on Eq. ( |5.10| ) . The two vertical lines 
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FIG. 3. The magnetogenesis region in the case 5 = 2 (i.e. increasing gauge coupling) and 
cp ~ 1. The exclusion plot is then given in terms of e and £. 



illustrate the conditions of Eqs. (|3.14j )- (|3.15|) . As in Fig. |] the two dashed lines correspond 
to the constraints coming from the inhomogeneous modes and from the condition e > £. 

The requirement that the spectra decrease at large distance scales implies, in this case 
that 5 < 2. The condition on the growth of the correlation function obtained in Eqs. Q4.ll )- 
( |4.12j ) imply 5 > 1/3. Thus the interesting physical range of Eq. ( [4.10| ) and ( |5.10| ) will be 
1/3 < S < 2. 

In Fig. |] the parameter space is illustrated for fixed values of £, i.e. £ ~ 1CT 11 . As in 
Fig. |3| the full and dashed curves correspond to the dynamo requirements imposed on Eq. 
( |5.10| ). The shaded area selects the allowed region in the space of the parameters where 
magnetogenesis is possible for 10 -11 < e < 10~ 6 . As in the case of Fig. ||] there are regions in 
the shaded area where values much larger than the magnetogenesis requirement are possible 
(see the dashed line in Fig. d). 

As it has been pointed out in deriving the theoretical bounds on the scenario the require- 
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FIG. 4. The magnetogenesis region is illustrated in the case of increasing gauge coupling in the 
(e, 5) plane. Notice that in the present example £ ~ 1CT 11 (for compatibility with the EW epoch) 
and ip ~ 1. 

ment £ > 10 -11 may be too restrictive since it excludes the variation of the gauge coupling 
at the electroweak time. To relax this assumption is possible and it would require a precise 
analysis of the dynamics of the EWPT in the presence of time varying gauge coupling. At 
the moment this kind of analysis is not available. Summarizing this illustrative discussion, 
there are regions in the parameter space of the models wlier all the theoretical constraints 
are satisfied and where magnetogenesis is possible. In particular, it is possible that the gauge 
coupling freezes prior to the electroweak epoch, leading still to magnetogenesis. 

VI. CONCLUDING REMARKS 

There are no reasons why the gauge couplings should be constant throughout all the 
history of the Universe. If they are allowed to change prior to the formation of the light 
elements they can lead to computable differences in the cosmological evolution. 

In the present paper the interplay between inflationary magnetogenesis and the evolution 
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of the Abelian gauge coupling has been addressed. In a phenomenologically reasonable model 
of inflationary and post- inflationary evolution the relaxation of the gauge coupling leads to a 
growth in the correlation function of magnetic inhomogeneities. Large scale magnetic fields 
are then generated. The evolution of the gauge coupling is driven, in the present context, 
by a massive scalar. 

Since the gauge coupling evolves significant changes in the evolution of the plasma can be 
envisaged. In the present investigation the ordinary MHD equations have been generalized to 
the case of time evolving gauge coupling but other effects could be envisaged. In particular, 
the generlization of the full kinetic approach to the case of time evolving "electron" charge 
would be of related interest. 

The value of large scale magnetic fields produced with this mechanism has been esti- 
mated. For a broad range of parameters the obtained values of the magnetic fields are much 
larger than the dynamo requirements. 

In the present investigation the main assumption has been that the only gauge coupling 
free to evolve is the Abelian one. Furthermore, the parameters of the model have been chosen 
in such a way that the gauge coupling is not dynamical by the onset of the electroweak phase 
transition. It would be interesting to relax both assumptions since they may lead to potential 
differences with the standard scenarios. 

Therefore, in many respects, the present investigation is not conclusive. At the same time 
it shows that acceptable models for the evolution of the gauge couplings can be obtained in 
a standard cosmological framework. 
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